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Abstract. We study (inverse) photoemission from systems with strong coupling of doped carriers to 
phonons. Using an adiabatic approximation, we develop a method for calculating spectra. This method is 
particularly simple for systems where the electron-phonon coupling can be neglected in the initial state, 
e.g., the undoped t-J model. The theory then naturally explains why the electron-phonon coupling just 
leads to a broadening of spectra calculated without electron-phonon coupling, without changing the dis- 
persion. This is in agreement with recent angle-resolved photoemission spectroscopy (ARPES) on undoped 
cuprates, and it supports the interpretation in terms of strong electron-phonon interaction. The theory 
also shows that for systems with strong electron-phonon coupling in the initial state, the result cannot in 
general be related to the spectrum obtained without electron-phonon coupling. 



PACS. 71.38.-k Polarons and electron-phonon interactions 
spectra - 71.10.Fd Lattice fermion models 



79.60.-i Photoemission and photoelectron 



1 Introduction 

Angle-resolved photoemission spectroscopy (ARPES) ex- 
periments have found evidence for strong electron-phonon 
interaction (EPI) and polaron physics in many materials 
like quasi-one-dimensional conductors [1,2], the mangan- 
ites [3], or the undoped high-T c cuprates [4-7]. The spec- 
tra show an incoherent broad feature whereas the quasi- 
particle peak at lower binding energy is strongly sup- 
pressed. In the case of the undoped high-T c cuprates, the 
dispersion of that broad peak matches the quasi-particle 
dispersion calculated in purely electronic models [6]. For 
an undoped t-J model with coupling of doped holes to 
optical phonons, numerical calculations of the ARPES 
spectra indeed showed broad features tracing the disper- 
sion of the quasi-particles in the original t-J model [8]. 
Similar observations in the manganites [3] and quasi-one- 
dimensional conductors [1,2] have been interpreted in anal- 
ogy with a single electron coupled to harmonic oscillators 
[9] and a related sum-rule for the first spectral moment. 
With respect to the manganites also the picture of the 
photohole seeing a frozen lattice has been used [10]. 

In this paper we address (inverse) photoemission spec- 
tra from systems with strong coupling of doped carriers to 
phonons. We develop a theory based on the adiabatic ap- 
proximation. This theory takes a particularly simple form 
if the EPI can be neglected in the initial state. This is 
the case for the undoped t-J model and the empty or full 
Holstein model. The spectrum can then be related to the 



spectrum of a model without EPI, and the effect of the 
EPI is essentially to broaden the spectrum. For systems 
where the EPI is important also in the initial state, the 
theory takes a more complicated form. It can then be re- 
lated to an average of spectra for distorted lattices without 
EPI. 

After giving general arguments in section 2 we discuss 
numerical results for the Holstein model (section 3) and 
the t-J model with phonons (section 4) as illustrating ex- 
amples. 



2 General considerations 

We consider a system that is modeled by the following 
Hamiltonian: 



H — H e i + H p h + H e 



(1) 



H e i (Hph) describes the purely electronic (phononic) part 
of the model whereas the interaction between electrons 
and phonons is given by H ep . 

The phonons are assumed to be harmonic in the ab- 
sence of EPI and the system is taken to be translationally 
invariant, so that we can write: 



H 



ph 
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Here, an individual phonon mode with frequency w q .„ has 
wavevector q and belongs to branch v. Its generalized co- 
ordinate and momentum are denoted by Q^.v and H^v 
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The EPI couples electronic degrees of freedom to the 
phonon coordinates <3 q ,„. We assume that this interaction 
vanishes for a certain electronic filling of the system which 
we will refer to as undoped in the following. The com- 
pletely empty or completely filled Holstein or Holstein-t- J 
model are examples for such undoped systems, see sections 
3 and 4. 

In addition, we assume that the phonon frequencies are 
small compared to the electronic energy scales defined by 
H e i. This justifies an adiabatic approximation [11] and we 
can first consider only H e i + H ep treating the phonon co- 
ordinates Qq M in H ep as c-numbers, i.e. as instantaneous 
parameters for the electronic problem. We denote the cor- 
responding eigenstates and eigenvalues by |^^ e (Q)} and 
Em" (Q) which are labeled by the number of electrons N e 
and other quantum numbers m. We use Q as a shorthand 
notation for the set of phonon coordinates. The phonon 
eigenfunctions are then obtained by solving 

(3) 

a Schrodinger equation with the effective potential 

C (Q) = K e (Q) + E 1<vQwQ-*,»- ( 4 ) 

The eigenenergies are e^ n where n stands for suitable 
phonon quantum numbers. In our approximation the eigen- 
states of if are (Q\s% n ) = 4> m " n (Q) \E%>(Q)), i.e. we as- 
sume that the electronic states |^^ e (Q)) do not mix. 

The ground-state in a system with N e electrons is then 
given by (Q|e£o) = 0oo (Q) \ E o"(Q)) with eigenenergy 
Sqq . In case of an undoped system with N° electrons there 
is no EPI and the phonon wavefunction just corresponds 
to the ground-state of H p h m . 

000° (Q) = nK,^) 1/4ex P(-^<?L/2). (5) 

The electronic ground-state and its eigenenergy are then 

independent of Q, and s q = E ° + J2q v w q,f /2- 

We now do (inverse) photoemission at zero tempera- 
ture by destroying (creating) an electron with momentum 
k and spin a in the ground-state of a system with N e 
electrons. Within the adiabatic approximation this can be 
described by considering the following Green's function 

<; ,T w = (<& iv> f — n .m &) (6) 

z~ {H-Soo ) 
= jdQ |rfQV V o*(Q)(< e (Q)IV' t x 

x(QI — * ^. |Q / )V'K e (Q0)^(Q / ) 

z l* 1 £ ao ) 

where ip = c^ a and J dQ = n q ,„ / dQ^ v . 

We proceed in analogy to Ref. [12] and neglect the ki- 
netic energy of the phonons in the resolvent in Eq. (6). 



Now H is diagonal in the phonon coordinates Q and one 
half of the integrations in Eq. (6) can be eliminated. This 
leads to the following approximation for the Green's func- 
tion [12]: 

G^ T (z) = JdQ |^'(Q)| 2 ff& ,=F (*,Q) (7) 

where 

z-{H(Q)-e 00 e ) 

(8) 

Sqq is the initial state's energy without any zero-point en- 
ergy contributions arising from the phonon kinetic energy. 
Finally, the corresponding spectral function is given by 

= JdQ\^(Q)\ 2 P^(^Q) (9) 

where 

pt T (^Q) = E K^ Tl (Q)IV>l< e (Q)>l 2 x 

m 

x<5(o;-(C eTl (Q)-4 e )) (10) 

after expanding ip\EQ e (Q)) in the adiabatic electronic ba- 
sis states |^ Tl (Q)). 

Equations (9-10) will turn out to be the key formula for 
interpreting ARPES spectra of undoped systems. To see 
this we observe that p^' T (w, Q) is the spectral function 
of the system without EPI for a given lattice distortion Q. 
If we assume that V N ' (Q) has a non-degencrate absolute 
minimum at Q m in the corresponding ground-state phonon 
wave- function will be localized around this point in coordi- 
nate space. If we approximate |0ocT (Q)| 2 ~ <KQ — Qmin) 
we find that the spectrum corresponds to the spectrum 
one obtains for the system with a frozen distortion Q m in 
in which there is no EPI. Analogously, in case of more 
than one (quasi-)degenerate minima of V^ f "(Q) we have 
to take the (weighted) superposition of the spectra cor- 
responding to the respective distortions. If we take into 
account the finite width of I^Jq (Q)| 2 it follows from Eq. 
(9) that the spectral features are broadened due to the 
Q-dependence of V r ^^ 1 (Q). We will consider a specific 
example in section 3. 

This analysis leads to our main conclusion. For the 
undoped system it follows from Eq. (5) that Q m in = 0. 
Consequently, the spectrum is just the broadened spec- 
trum of the same system without EPI (H ep = 0). The 
dispersion of the (H ep = 0)-quasi-particle peak shows up 
in the k-dependence of the broadened peak in the low 
binding energy part of the spectra. This approach is par- 
ticularly useful for the undoped system as it allows state- 
ments about the spectrum for H ep ^ from the knowledge 
of the spectrum for H ep = 0. For the doped system Q m i n 
is typically non-zero if the EPI is strong. The spectrum 
/ , k ? CT =F ( tJ ' Qmin) can then be very different from the spec- 
trum for Q = 0. Therefore, even if the (H ep = 0)-spectrum 
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is known, in general no information about the spectrum 
of the system with EPI can be deduced. In this context 
we notice that there may be several degenerate minima at 
Qmm 7^ such that the ground-state has an undistorted 
lattice in the sense of a vanishing expectation value of Q. 
The spectrum for H ep ^ 0, nevertheless, corresponds to 
a superposition of the spectra for Q-values around the 
minima Q 

////// 

The numerical calculation of spectral functions using 
the approximation in Eq. (9) can be quite efficient com- 
pared to other methods typically used to obtain spec- 
tra. In general, it is much easier to obtain the spectrum 
Q) for a system without EPI but a given distor- 
tion Q than for a system with EPI. If the initial state's 
phonon wavefunction ^^(Q) is known (e.g. in the case 
of an undoped system) its square can be used as a weight 
function in a Monte Carlo integration over the phonon co- 
ordinates in Eq. (9). The numerical effort is independent 
of the strength of the EPI. On the other hand, if e.g. exact 
diagonalization is used to obtain directly the ARPES spec- 
trum of a system with EPI, calculations become computa- 
tionally very demanding with increasing coupling strength 
as the truncated phonon Hilbert space grows larger and 
larger. 

In order to improve the present approximation we have 
to include also the kinetic energy of the phonons in the 
resolvent of the Hamiltonian. Then, iP\Eq s (Q)) in Eq. (6) 
must be expanded with respect to both electronic and 
phononic basis functions in the adiabatic approximation. 
From this, one obtains the following expression for the 
spectral function: 



<y =F M=E <y ( w -& =F1 - e ~)) 



(ii) 



Each eigenstate of H in the (A e =F l)-electron sector repre- 
sents a possible final state. It contributes to the spectrum 
at its eigenenergy (shifted by the ground-state energy of 
the system with N e electrons). The intensity is propor- 
tional to the squared overlap of final and initial state, i.e. 
the ground-state of the system with N e electrons plus an 
additional hole or electron. Two conditions must be ful- 
filled for this overlap to be large. 1 

i) As the initial phonon wavefunction <^>oo e (Q) is local- 
ized around the minimum Q m in of V^^Q) and has no 
nodes, the final phonon wavefunction 4>mn Tl (Q) must have 
a large and slowly varying amplitude in this region, too. 
This will be the case for final states with energies e^f 1 w 
Kra 6 ^ '(Qmin) • For smaller energies the region around Q m in 
is classically forbidden and the amplitude of the final pho- 
non wavefunction becomes exponentially suppressed, where- 
as for larger energies the kinetic energy increases and the 
wavefunction oscillates faster. In both cases the integrated 



overlap of initial and final phonon wavefunction becomes 
smaller again. 

ii) The electronic matrix element (E^ 1 (Q)\tp\E^ (Q)) 
must be large around Q m i n . It is sufficient to consider its 
value only in this region as the initial phonon wavefunc- 
tion is small elsewhere. 

Altogether this leads to the following picture: For a 
system without EPI but with a given lattice distortion 
Qmm the spectrum consists of (5-functions at the energies 
Em" Tl (Qmin) ■ If the EPI is switched on, spectral features 
with large intensities will still appear at similar energies 
and with similar relative weight but they will be broad- 
ened by phonon sidebands. The quasiparticle's dispersion 
and weight, however, can be strongly altered by the EPI. 
In general, the effective phonon potential V^ VeTl (Q) cor- 
responding to the electronic ground-state in the system 
with N e =p 1 electrons has minima at Q ^ Qmin- The 
ground-state phonon wavefunction 0^[) e=Fl (Q) is localized 
around these minima. Consequently, there is only little 
overlap with the phonon wavefunction in the initial state 
which peaks around Q m i„ and in the spectrum the peak 
lowest in binding energy has only very small weight. 

It is interesting to discuss the problem above in terms 
of a sum-rule concerning the first moment (center of grav- 
ity) of the spectrum. For the undoped system with N® 
electrons that has no phonons excited in the initial state, 
one can show under rather general assumptions that the 
first moment of the (inverse) photoemission spectrum does 
not depend on the strength of the EPI. Let us consider 

the case when A k y (ui) has only one peak for a given 
k in the absence of EPI, e.g. if there is only one band 
and no electron-electron interaction. Turning on the EPI 
then broadens the peak, but due to the sum-rule the in- 
teraction is normally not expected to drastically change 
the spectrum in other respects. If the electrons interact 
with themselves, however, already for systems without 

EPI A k " a (uj) usually has several peaks for a given k 
and the first moment does not correspond to the quasi- 
particle energy. Then, the sum-rule is not able to tell us 
how the different peaks are broadened and cannot be used 
to argue for prominent features in the spectrum dispers- 
ing approximately like the quasi-particles in the system 
without phonons. 



3 Holstein model 

In the following we demonstrate the validity of the argu- 
ments given above with several examples. First, we con- 
sider the one-dimensional A-site Holstein model with pe- 
riodic boundary conditions for which the electronic part 
of the Hamiltonian just contains nearest-neighbor hopping 
with amplitude t: 



H el = -2t^cos(fc) c{ a c k a , 



(12) 



1 In the following discussion we assume a non-degenerate 
minimum of 1q JVc (Q) for simplicity. The arguments can be eas- 
ily generalized to the case of (quasi-)degenerate minima. 



where c\. creates an electron with momentum k and spin 
a. The electrons are coupled to dispersionless phonons 
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with frequency co p h as described by 



g_ 

N 



£ 

q,k,(7 



ph 



3 C k+q,a C k,a- 



(13) 



The strength of the interaction is given by g. The (q = 0)- 
phonon mode couples to the total number of electrons, 
N e . We can therefore write H = Hi + H 2 , where 



H, 



9 



i 



V^Qo + -z (il 2 + Lu 2 ph Q 2 a ) 



(14) 



can be solved exactly [9]. The spectral function A^" a T ^ (lu) 
obtained for Hi only needs to be convoluted by the known 
result for Hi, 

oo i 

A N "*W (w) ^-^^( W + (1 T 2N e )aw ph - lu> ph ), 



/=o 



where a = (g/ui p h) 2 /N , in order to get A^° a T (uj) 



(15) 
as de- 



fined in Eq. (11). We will therefore restrict our discussion 



to A 



JVe,=F(2) 
k,<j 



(u) in the following. 



We specifically consider a two-site system (N = 2) 
[13] and calculate the inverse photocmission spectrum for 
creating an electron in both the empty (undoped) system 
and the system that already contains an electron of oppo- 
site spin. There is then only one phonon coordinate Q v in 
Hi which we treat as a parameter in solving the part of 
Hi coming from H e i + H ep . In the one-electron sector one 
obtains the two eigenenergies 

= Ty/t*+u ph g*Ql, (16) 

whereas there are four eigenergies in case of two electrons 
with opposite spin: 



^0/3 



= T 2^t 2 +Lo ph g 2 Q 2 7n </ 2 =2 =0. (17) 



For numerical calculations we choose t = 1, ui p h = 0.1, and 
g = 0.6. As Lu p h/t is small our adiabatic approximation is 
justified. The effective potentials one obtains by adding 
the harmonic potential u> ph Q 2 /2 to the eigenenergies in 
Eqs. (16) and (17) are shown in Figs. 1 and 2, respectively. 

We first consider the approximation when the kinetic 
energy of the phonons is neglected in the resolvent of the 
Hamiltonian and for which the spectral function is given 
by Eq. (9). The neglected terms are proportional to u> p h so 
we cannot expect to resolve fine-structure in the spectra on 
that order. But it turns out that this approximation still 
describes the overall broadening correctly on a larger scale 
proportional to ^/uJpTi [12]. The phonon wavefunction in 
the initial state is known exactly for the undoped system 
, see Eq. (5). Here 



9 oo 



(Q*) = K^A) (1/4) cM-u P hQl/2). 



(18) 



In the system with one electron the lowest effective po- 



tential Vq e 1 (Q ir ) has two minima at Q± around which 
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Fig. 1. Effective potentials in the two-site Holstein model with 
one electron as functions of the phonon coordinate Q n for i = 1, 
LOph = 0.1, g = 0.6. Some selected phonon wavefunctions are 
also shown, see text. 
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Fig. 2. Effective potentials in the two-site Holstein model with 
two electrons with opposite spin as functions of the phonon 
coordinate for t = 1, u) p h = 0.1, g = 0.6 are shown in 
the left panel together with the phonon wavefunction ifioo^ 1 
(dotted line). Right panel: A^^ +{2) (cu) rotated by 90°. 



N„ 



we can treat it as harmonic potential with renormalizcd 
phonon frequency ZJ = ^J(d 2 V~Q 
to the approximation 



1 /3QI)\q . This leads 



h N c = l 



v l=± 



1/4 



exp(-C7(Q 7r - Q t ) 2 /2). 



(19) 

Following Ref. [12], we can expand the argument of the 
S- functions in Eq. (9) up to first order in Q n around Q ± . 
If we further assume that the electronic matrix elements 
vary only weakly around Q± the integration over in 
Eq. (9) can be eliminated and we obtain the following 
result: 



2 ^\El) 



-(a 



o 2 



(20) 
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Fig. 3. Spectrum for adding an electron to a two-site Holstein 
model with one electron of opposite spin (t = 1, uj p h = 0.1, 
g = 0.6). The approximation A^q^ + (w) (Eq. (20)) is shown 
as dashed line together with the spectrum obtained from exact 
diagonalization (thin line). Both spectra have been convoluted 
with a Lorentzian (FWHM=0.1). 
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-2-1 1 2 

CO 

Fig. 4. Spectra for adding an electron to an empty two-site 
Holstein model (t — 1, uj p h = 0.1, g = 0.6). The approximation 

~ iV^— + (2) 

a ( w ) (^1- (9)) thin lines) and results from exact diag- 

onalization (bold lines) are shown for both k — (solid lines) 
and k = -K (dashed lines). All spectra include a Lorentzian 
broadening (FHWM=0.06 eV). 



where a m = V^\q ± ~ e x m and b m = (dV^/dQ^)\^ ± . This 
approximation is shown in Fig. 3 for k = together with 
the spectrum obtained from exact diagonalization. The 
agreement with the Gaussian lineshape predicted by Eq. 
(20) is very good. Only to resolve the fine-structure on a 
scale given by the phonon frequency uj p h one has to go 
beyond the present approximation. As indicated schemat- 
ically by the arrows in Fig. 2 the spectrum can indeed be 
understood as the spectrum of electrons in a system with- 
out EPI but a given distortion Q±. The broadening is due 
to the finite width of the phonon wavefunction in the ini- 
tial state. There is no structure in the spectrum arising 
from the highest effective potential as the corresponding 
electronic matrix element vanishes. 

In the case of inverse photoemission from the empty 
system the phonon wavefunction in the initial state is lo- 
calized around Q n = 0. The slope of both effective po- 
tentials V^.q \(Qir) vanishes at this point. An approxi- 
mation analogous to Eq. (20) would therefore result in 

A^ e ~°' + ^ 2 \u)=p^ e ~°' + {w, Qir = 0), i.e. the unbroadened 
spectrum of the undistorted system without EPI. If we 
evaluate Eq. (9) without any further approximation we 
obtain the spectra shown in Fig. 4. The comparison with 
results from exact diagonalization shows that this approx- 
imation cannot reproduce the fine-structure and for k = 
does not give spectral weight above uo = — 1 (except from 
broadening), but it still gives the right order for the broad- 
ening of the peaks. 

In order to understand also details of the spectra we 
have to use Eq. (11) which results from making no ap- 
proximation other than the initial adiabatic one. One has 
to solve Eq. (3) for each effective potential to obtain the 
phonon eigenfunctions needed in Eq. (11). The result- 
ing inverse photoemission spectra for creating an electron 
with momentum k = or k = it in the empty system 
are shown in Fig. 5. They are practically indistinguish- 



100 



0.0001 




0.01 



0.001 r 



N — +(2) 

Fig. 5. Spectra A k " a ' (u>) for adding an electron to an 
empty two-site Holstein model (t = 1, tu ph — 0.1, g = 0.6) 
as given by Eq. (11). The spectra for both k — (solid line) 
and k = it (dashed line) are shown with Lorentzian broadening 
(FWHM=0.01). Observe the logarithmic intensity scale. 



able from those we obtained using exact diagonalization. 
This shows that the adiabatic approximation works very 
well for the chosen parameters. For the electronic matrix 
element in Eq. (11) one needs the Q T -dependent electronic 
eigenstates corresponding to the eigenenergies in Eq. (16): 



(21) 

with 



(sign(Q T )) Q 



y/(t+ (-1)"V* 2 + u ph g 2 Ql) 2 + u ph g 2 Ql 



, (22) 



where |0) and \ir) denote the k = and k = it one-electron 
states and m = 0, 1. 
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Let us consider e.g. the case where we create an elec- 
tron with momentum k = in the empty system. The ini- 
tial electronic state is then |0). The initial state's phonon 

N o =0 

wavefunction OO ° (dashed line in Fig. f ) is a bell-shaped 
Gaussian centered around = 0. This effectively lim- 
its the integration over in Eq. (11) to a small region 
around the origin. Here \E*' =1 ) w |0) and |£f e=1 ) « \ir). 
Therefore, the relevant electronic matrix element in this 
example is (Eq" =1 \0). It has even parity with respect to 
Qn- So, according to Eq. (11) final states with the elec- 
tronic configuration |£^ c_1 ) and a corresponding even- 
parity phonon wavefunction </>^ =1 that strongly overlaps 
with (^^f -0 give rise to large spectral intensity. We show 
in Fig. 1 the even-parity phonon wavefunction (pQ^r 1 with 
the largest overlap offset along the ordinate by its eigen- 
ergy. It has a sizable and slowly varying amplitude around 
Qtt = because its eigenergy is close to the local value 
of the effective potential which at Q n = equals the 
(k = 0)-eigenenergy —t of the system without EPI. There- 
fore, the peak with largest weight appears around this en- 
ergy in the spectrum (solid line in Fig. 5). The sidepeaks 
arise from final states with even-parity phonon wavefunc- 
tions with lower or higher eigenenergies whose overlap 

with (t>oQ~° decreases. Figure 1 also shows the even parity 
ground-state phonon wavefunction ^^f -1 in the double- 

N°— 

well potential. Clearly, its overlap with c/> 00 c is very small 
leading to a strongly suppressed quasi-particle peak in 

Ak=oa' + ^- Because of the large dimensionless EPI con- 
stant A = g 2 /(ui p ht) = 3.6 we are well in the polaronic 
regime. 

If, on the other hand, an electron with momentum k = 
7r is created in the empty state one finds using similar 
arguments as before that final states leading to a large 
spectral intensity must have the electronic configuration 
\E^"~ 1 ). Their phonon wavefunction must strongly 

overlap with (j)^ ~° and be of even-parity. In this case the 
lowest energy phonon wavefunction 0^f =1 in the upper 
effective potential (shown in Fig. 1) has the largest overlap 
because its eigenenergy is closest to E^'~ 1 (Q 7r = 0) = +t. 

— I f 2~) 

Therefore, A k ^ Q a (u>) (dashed line in Fig. 5) shows a 
prominent peak at w ~ +t. 

As (£f e=1 |0) ((E^" =1 \n)) only vanishes completely at 
Qtt = the spectrum for k = (k = tt) also shows weak 
structures around u) = +t (u) = —t) where the coupling is 
now to phonon wavefunctions of odd parity. The density of 
coupling phonon states is different around ui = — t and ui = 
+t. E.g., the fact that the upper effective potential has a 
minimum around = results in an asymmetric shape 

N (> +(2) 

of A k y (ui) around u> = +t as no phonon cigenstates in 
this effective potential can have eigenenergies below +t. 

4 t-J model with phonons 

As a second example we study the one-dimensional TV-site 
Holstein-i- J model with periodic boundary conditions. This 
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Fig. 6. Spectral functions A k h a ' (w) for creating a hole in 
the undoped 4-site Holstein-t- J model (t = 1, J = 0.3, io ph — 
0.1, g = 0.8). k = 0: solid line, k = ±7r/2: dashed line, k — 
7r: dotted line (Lorentzian broadening: FWHM=0.01). Arrows 
show positions and weights of corresponding peaks for g — 0. 

model also includes electron-electron interactions. The elec- 
tronic part of H is given by the usual t-J Hamiltonian 

+ 7vj( s , Si+1 _^±i), 

i 

where cj a creates an electron with spin a on site i if this 
site was previously empty, n, = cJ i0 -C i o ., and Sj is a 
spin-i operator. Besides the hopping t there is also an ex- 
change coupling parameter J. As in the Holstcin-model we 
consider an interaction with dispersionless phonons where 
the coupling is now to empty sites (holes): 

H -p = J^Y1 V^Q^ 1 - n iV q °- ( 24 ) 

1,3 

The system with one electron per site corresponds to the 
undoped case where the EPI vanishes. The (q = 0)-phonon 
mode can be treated separately again, the only difference 
being that the coupling is now proportional to the total 
number of empty sites Nh, not to the total number of 
electrons N e . 

For numerical calculations we consider a 4-site system 
with t = 1, J = 0.3, ujph = O.f, and g = 0.8. The photoe- 
mission spectra for destroying an electron with momentum 
k and spin a were obtained using exact diagonalization 
with up to 200 phonons per basis state for solving the 
problem without the (q = 0)-mode and subsequent convo- 
lution with A N h=°-M(<jj) (Eq. (15) with N h = instead 
of N e ). The results are displayed in Fig. 6 together with ar- 
rows indicating the peaks in the corresponding spectra for 
g = 0. Without EPI there is only one peak both for k = 
(at —0.55) as well as for k = tt (at f .45). For k = ±tt/2 the 
spectrum has two peaks at —1.139 and 2.339 as a result 
of the electron-electron interaction. 
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Fig. 7. Inverse photoemission spectra (k = (3, 1)tt/5) for the 10-site t-J model with one interacting phonon mode (t = 1, 
J = 0.4, uj p h = 0.1, <7(7r,jr) = 1.6): (a) from the undoped system, (b) from the 10%-doped system. The corresponding spectra for 
systems without EPI are shown by dashed lines with their amplitude flipped for clarity. A Lorentzian broadening of FWHM=0.04 
has been applied. 



Again, the dispersion found in the system without EPI 
is traced quite accurately by a broad peak in the case of 
strong EPI. The spectrum for k = ±7r/2 also illustrates 
our comments at the end of section 2 on the sum-rule 
concerning the first spectral moment. Since the spectrum 
has two peaks for k — ±7r/2, the sum rule cannot tell us 
how the peaks are broadened. Many other spectra would 
also have been consistent with the sum-rule, e.g. spectra 
where the peaks are shifted. The arguments based on the 
adiabatic approximation, however, show that both peaks 
should be broadened with their individual center of gravity 
remaining roughly unchanged in agreement with the exact 
calculations. 

Finally, we consider the t-J model in two dimensions 
on a tilted 10-site cluster with periodic boundary condi- 
tions. To simplify calculations we assume that the EPI is 
described by Eq. (24) but with a q-dependent coupling 
constant <? q . In the following we choose ff( W)ir ) = 1.6 and 
5 q = for all other q ^ (tt, tt) so that effectively there 
is only one phonon mode that interacts with the elec- 
trons. The other parameters are t = 1, J = 0.4, and 
uj p h =0.1. Figure 7 shows the inverse photoemission spec- 
tra for k = (3, l)7r/5 from both the undoped (Fig. 7(a)) 
and the 10%-doped system (Fig. 7(b)). The EPI has been 
switched on and off (solid line vs. dashed line with flipped 
amplitude). 

The spectra from the undoped system confirm again 
our general expectations from section 2. The EPI basically 
broadens the structures in the original spectrum. This in- 
cludes the quasi-particle peak at low binding energies. In 
contrast, the spectrum from the doped system changes 
quite differently when the EPI is switched on. Although 
the spectrum again develops several broad features they 
cannot be related anymore in a simple way to the struc- 
tures in the spectrum found without EPI. We have also 
varied k and found that the dispersion of the broad fea- 
tures is different from the quasi-particle dispersion in the 
system without EPI. According to Eq. (9) the spectra 



rather correspond to broadened versions of spectra one 
would obtain in a purely electronic, but distorted system. 



5 Conclusions 

We have introduced an adiabatic approximation for calcu- 
lating ARPES spectra from systems with strong coupling 
of doped carriers to phonons. The effective phonon poten- 
tial for the initial state is calculated as a function of the 
phonon coordinates and its minima are found. We show 
that the spectrum with electron-phonon interaction (EPI) 
is then related to a broadened average of spectra with- 
out EPI calculated for distorted lattices corresponding to 
the minima of the effective potential. We also studied the 
additional approximation of neglecting the kinetic energy 
of the phonons in the resolvent of the Hamiltonian corre- 
sponding to the ARPES Green's function. The spectrum is 
then expressed as superposition of spectra from distorted 
lattices without EPI, using the square of the initial state's 
phonon wave function as a weight function (see Eq. (9)). 

In either form, the theory becomes particularly simple 
if the EPI can be neglected in the initial state. The phonon 
wave function is then centered around the undistorted lat- 
tice, and the spectrum with EPI can be directly related 
to the (broadened) spectrum without EPI. In the case of 
strong EPI in the initial state, the minima of the effective 
potential correspond to distorted lattices. The spectrum 
with EPI is then related to the (broadened) spectra of dis- 
torted lattices without EPI. Therefore, the knowledge of 
the spectrum without EPI for the undistorted lattice is in 
general not very informative with respect to the spectrum 
with EPI. 

Our results support the interpretation of ARPES on 
undoped high-T c cuprates [6,7] and explain why in nu- 
merical calculations [8] the quasi-particle dispersion from 
purely electronic models shows up almost unchanged in 
the dispersion of incoherent features in the spectra ob- 
tained with EPI. 
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